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摘要 : 本 文 简要 介绍 新 的 对 称 形式 的 质点 动力 学 方程 对 天 体 动 力学 理论 的 系统 性 改进 。 首 
先 ， 对 于 开放 的 多 体系 统 ， 无 法 找到 近似 程度 非常 高 的 惯性 系 , 在 传统 理论 中 不 得 不 引入 惯 
性 系 的 近似 ， 而 对 称 新 方程 因为 可 以 直接 适用 于 任意 的 平 动 参考 系 而 避免 了 惯性 系 的 近似 ， 
从 而 可 以 提高 理论 预言 的 精度 。 其 次 ， 对 于 束缚 的 多 体系 统 ， 传 统 理论 的 动力 学 应 用 是 先 引 
入 质心 参考 系 , 在 质心 参考 系 中 应 用 牛顿 第 二 定律 , 然后 通过 坐标 变换 再 转化 到 实体 参考 系 ， 
比如 太阳 系 的 行星 摄 动 方程 。 但 是 ， 应 用 对 称 新 方程 则 可 以 一 步 到 位 推导 得 到 行星 摄 动 方 
程 。 最 后 ， 如 果 进 一 步 考 虑 行星 受到 临时 推力 或 者 冲击 力 ， 甚 至 为 了 进一步 提高 计算 精度 进 
而 考虑 来 自 束 缚 系统 外 的 作用 力 , 则 一 个 可 以 登 加 非 微 扰 作 用 力 的 行星 摄 动 的 修正 方程 在 本 
文 得 到 了 确立 。 
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Abstract: This article presents a brief introduction to a Systematic improvement 
in celestial dynamics theory through the introduction of a new symmetrical form of 
particle dynamics equation. For open multi-body systems where it is challenging to 
find a highly accurate inertial reference frame, approximations of inertial frames 
need to be introduced. However, the new symmetric equation can be applied to any 
translational reference frame, thus avoiding the need for inertial reference frame 
approximations and enhancing the accuracy of theoretical predictions. In the case 
of bound multi-body systems, the traditional approach involves introducing a center 
of mass reference frame, applying Newton’s second law in this frame, and then 
transforming it back to the actual reference frame, such as the perturbation equation 


of the solar system s planets. Conversely, applying the new symmetric equation 


allows for the direct derivation of planetary perturbation equations. Furthermore, 
by considering temporary thrust or impact forces acting on planets, or to further 
enhance computational precision, even considering any external forces acting on the 
bound system, then a new correction equation has been established for the planetary 
perturbation equation that can be further imposed with non-perturbative 
interactions. 

Keywords: Celestial dynamics Non-perturbative correction New symmetric 


equation Planetary perturbation equation 


* 本 文系 中 科 院 知识 创新 工程 基金 项 目 (项 目 编 号 ;: KJCX2-YW-W10) 的 研究 成 果 之 一 。 


1 引言 


传统 天 体 动力 学 的 理论 基础 是 牛顿 第 二 定律 , 而 牛顿 第 二 定律 的 适用 条 件 是 
惯性 参考 系 。 仅 以 此 为 基础 ， 实际 上 可 以 严格 推导 出 一 个 对 称 形式 的 质点 动力 学 


方程 “在 理论 上 假定 宇宙 中 存在 某 一 个 惯性 系 >, 待考 察 的 运动 天 体 为 p， 任 


意 选 取 一 个 实际 的 天 体 为 参考 物 o。, 原则 上 我 们 可 以 运用 牛顿 第 二 定律 求 出 这 两 
个 天 体 相对 于 惯性 系 的 加 速度 ， 分 别 为 
F|,=m,al, s 
F| =m, ,a| > 
为 了 展现 形式 上 的 对 称 性 ， 上 和 式 又 可 以 写成 ， 


(1) 


(2) 


a 根据 牵连 加 速度 、 相 对 加 速度 、 绝 对 加 速度 之 间 的 人 
利 略 变换 关系 ”， 

al,s 二 人 a (3) 
现在 以 o 为 参考 原点 ， 建立 相对 宇宙 空 人 
标记 为 O (为 了 区 分 , 这 里 约定 平 动 参考 系 及 其 参考 质点 分 别 由 对 应 的 大 写 英 文 
字母 和 小 写 英 文字 母 标 记 ) 。 考虑 到 惯性 0 间 背 景 无 旋 
转 的 ， 则 有 


(al,,). =a|,, (4) 
我 们 将 式 (2) (3) (4) 进行 联 立 ， 得 
， 攻 (5) 
m m 0 


此 即 新 近 引 入 的 对 称 形式 的 质点 动力 学 方程 (下 面 简称 为 对 称 新 方程 》， 其 中 强 
含 的 物理 意义 和 多 重 推导 已 经 有 详细 的 讨论 ”” 。 本 文 则 重点 展示 其 在 天 体 动 力 
学 的 应 用 优势 。 显 然 ， 对 于 开放 的 多 体系 统 ， 应 用 午 顿 第 二 定律 的 传统 策略 就 是 
把 宇宙 天 体 中 的 茶 一 个 参考 系 近似 为 惯性 系 。 比 如 研究 太阳 和 月 球 的 相对 运动 时 
把 日 心 参考 系 近似 为 惯性 系 ”“， 这 种 方案 ， 即 便 近 似 程度 再 高 也 注定 会 带 来 一 定 
的 误差 。 针 对 这 种 由 于 理 双 本 身 的 适用 范围 带 来 的 误差， 应 用 对 称 新 方程 来 直接 
求解 非 惯性 系 下 的 天 体 动力 学 ， 既 简化 了 步骤 ， 又 取消 了 惯性 系 近 似 带 来 的 误 
差 ， 从 而 精确 了 计算 结果 ”“。 当然 在 传统 理论 中 ， 对 于 束缚 的 系统 ( 即 忽略 外 
部 受 力 的 近 孤 立 系 统 ) ， 还 存在 一 种 更 精确 求解 茶 个 天 体 动力 学 的 方案 。 比如 
考察 整个 太阳 系 内 某 一 个 行星 受到 的 除 太 阳 之 外 的 其 他 行星 的 引力 扰动 , 存在 闭 
名 的 行星 摄 动 的 微分 方程 。 接 下 来 本 文 即 利用 对 称 新 方程 重新 考察 行星 摄 动 的 
微分 方程 的 推导 和 修正 。 


2 传统 行星 摄 动 方程 的 验证 性 推导 


以 近 孤 并 的 任意 多 体系 统 为 例 , 应 用 对 称 新 方程 重新 推 时 行星 摄 动 微分 方程 
“。 在 经 典 的 太阳 系 n 体 问 题 中 , 设 n 为 太阳 ， 其 他 n-1 个 天 体 为 行星 。 直 接 引 


入 太阳 实体 参考 系 , 直接 (i<n) 在 日 心 参考 系 中 的 坐标 为 (%,y,z)。 
阁 第 ;个 


根据 对 称 新 方程 (5)， 任 意 第 i 个 行星 相对 日 心 参 考 系 (参考 质点 为 4n， 故 对 应 
平 动 参考 系 为 N ) 的 加 速度 满足 


设 第 i 星 


个 行星 


(6) 


把 这 个 矢量 方程 分 解 成 (x,y,z) 三 个 方向 的 分 量 方程 ,我 们 取 x 方 向 的 方程 进行 接 


下 来 的 推导 
FI F 
m, m 


本 节 和 暂且 假定 只 考虑 太阳 系 内 引力 相互 作用 , 则 根据 万 有 引力 定律 , 我 们 可 以 得 


到 灭 | ,FF ， 即 
i#j Gmm,(x,—x, 
Fl= > el 
j=1...n 六 
(8) 
mlG (2 一 
| 二 ， 本 
J=1 万 
只 要 将 具体 受 力 公式 (8) 代入 对 称 新 方程 的 分 量 形式 (7) 式 ， 立 刻 就 得 
Gm 一) 四 G 和 一 为 四 Gm (Xj—% Gm, ( wl Gm, 
eR Re 
j=l.n A; j=l A, j=1.n-l Ai j=l Ai 
站 
Gm (Xj—% m, +m GMX, 
- 用 Ai r | 
按照 惯例 ,现在 引入 直接 包含 第 三 方 天 体 参数 (所 有 带 指标 j 的 参数 ) 的 摄 动 函 
数 RR594， 其 定义 为 


十 yy 十 ZZ， 
neon [Le 10) 
Ai 广 
将 (10) 代入 (9) 式 可 得 
(m, +m, )x i OR, 
E+G— =  》 一 小 (11) 
r J=1...n—l Ox; 


i 


对 y 和 z 方 辐 的 方程 也 可 以 进行 和 上 式 类 似 的 推导 。 此 即 为 天 体力 学 中 行星 摄 动 


的 基本 微分 方程 ”。 因 此 ， 对 称 新 方程 (5) 在 经 典 力学 框架 下 的 正确 性 由 此 再 次 
得 到 了 验证 ,但 是 由 于 对 称 新 方程 对 惯性 系 依赖 的 规避 , 使 得 推导 过 程 一 步 到 位 ， 
具体 解 题 过 程 可 以 得 到 有 效 的 简化 。 


3 行星 摄 动 方 程 的 非 微 扰 修正 


在 上 一 节 的 推导 中 , 实际 上 假定 n 个 天 体 都 只 受到 系统 内 其 他 天 体 的 引力 相 
互 作用 ， 但 是 在 实际 情况 下 ， 天 体 可 能 受到 其 他 非 引 力 相互 作用 的 影响 。 更 需 指 
出 的 一 点 是 , 行星 摄 动 方程 的 推导 经 由 质心 参考 系 来 引入 牛顿 第 二 定律 ,看 似 精 
确 ， 但 是 实质 上 是 在 更 高 层次 上 对 惯性 系 作 了 近似 ， 即 把 所 考察 的 整个 天 体 束 
缚 系统 的 质心 参考 系 近似 为 了 惯性 系 ， 因 此 ， 在 特殊 情况 下 为 了 更 高 的 精度 ， 有 
必要 考虑 来 自 n 体 束缚 系统 之 外 的 作用 力 。 有 了 对 称 新 方程 (5) ， 这 个 问题 可 以 
得 到 简洁 严格 的 理论 解决 方案 。 

假定 第 i 个 天 体 , 受到 的 引力 相互 作用 标记 为 | ， 而 受到 的 任意 来 源 的 非 


引力 相互 作用 标记 为 fa|， 受到 的 束缚 系统 之 外 的 作用 力 统一 标记 为 fa4|， 
此 ， 第 (6) 式 可 以 改写 为 
(fal+ faal+ fenl,) (fal, + fal, + fa 


a (12) 


m, m, 
重复 从 方程 (7) 至 方程 (11) 的 推导 ， 不 难得 到 x 方向 的 动力 学 满足 
$d Sons. 是 oa pad 0 


3 3 3 
| A; 天 ET 


同样 引入 摄 动 函数 Rj 之 后 ， 登 加 了 非 微 扰 作用 力 之 后 的 行星 摄 动 微分 方程 修正 
为 


区 +G 


(m+m, )X > 和 A 的 


rn | Ox, m, m, 
其 中 对 于 非 引 力 相 互 作用 ， 可 以 是 火箭 的 推动 力 ， 也 可 以 是 爆炸 等 破坏 力 ， 
此 ， 上 述 方程 (14) 在 人 为 改变 小 行星 轨道 的 理论 计算 中 可 能 有 一 定 的 作用 。 


4 总 结 与 展望 


本 文 调用 了 新 近 引 入 的 对 称 形 式 的 质点 动力 学 方程 研究 了 其 在 天 体 动 力学 
的 理论 改进 ,对 称 新 方程 可 以 根据 理论 上 存在 的 惯性 系 从 牛顿 第 二 定律 严格 地 推 
导 得 到 , 本 文 证 明 其 在 多 体系 统 的 应 用 方面 具备 不 可 忽视 的 优势 。 对 于 一 般 的 开 
放 性 多 体系 统 , 可 以 因为 避免 惯性 系 的 近似 而 提高 计算 的 精度 。 对 于 束缚 的 近 孤 
并 多 体系 统 ， 对 称 新 方程 可 以 有 效 简化 对 行星 摄 动 方程 的 推导 ; 在 此 基础 上 ， 
进一步 考虑 行星 的 临时 新 增 受 力 或 者 深究 来 自 系 外 的 作用 力 ， 则 得 到 了 行星 摄 
动 的 非 微 扰 修 正方 程 。 
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